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ABSTRACT 

Numerical-relativity simulations indicate that the black hole produced in a binary merger can recoil 
with a velocity up to u max ~ 4,000 km/s with respect to the center of mass of the initial binary. This 
challenges the paradigm that most galaxies form through hierarchical mergers, yet retain supermassive 
black holes at their centers despite having escape velocities much less than u max . Interaction with a 
circumbinary disk can align the binary black hole spins with their orbital angular momentum, reducing 
the recoil velocity of the final black hole produced in the subsequent merger. However, the effectiveness 
of this alignment depends on highly uncertain accretion flows near the binary black holes. In this Letter, 
we show that if the spin Si of the more massive binary black hole is even partially aligned with the 
orbital angular momentum L, relativistic spin precession on sub-parsec scales can align the binary black 
hole spins with each other. This alignment significantly reduces the recoil velocity even in the absence of 
gas. For example, if the angle between Si and L at large separations is 10 degrees while the second spin 
S2 is isotropically distributed, the spin alignment discussed in this paper reduces the median recoil from 
864 km/s to 273 km/s for maximally spinning black holes with a mass ratio of 9/11. This reduction will 
greatly increase the fraction of galaxies retaining their supermassive black holes. 

Subject headings: black hole physics — gravitational waves — galaxies: evolution 



1. INTRODUCTION 

Observations suggest that most galaxies host super- 
massive black holes (SBHs) at their centers whose 
masses are ti ghtly correlated with properties of their 
host s pheroids (IMagorrian et al.lll998HFerrarese fc MerrittJ 
2000| iTremaine et al.ll2002|) . If galaxies form through hi- 
erarchical mergers, their SBHs may form from the merger 
of the smaller SBHs in their progenitor galaxies. The final 
stage of these black-hole mergers involves highly curved, 
dynamical spacetime that can only be simulated with fully 
nonlinear numerical r elativity (NR). Following a major 
breakthrough in 2005 dPretorius 1 120051 : ICampanelli et al.l 
120061 iBaker et al.l ImaT numerical relativists can now 
accurately determine the anisotropic emission of gravita- 
tional waves during the final stage of black-hole mergers. 
When gravitational waves are preferentially emitted in one 
direction during a merger, conservation of linear momen- 
tum requires that the final black hole produced in that 
merger recoil in the opposite direction. These recoil ve- 
locities or "kicks" can approach 4,000 km/s for maximally 
spinni ng mergers (jCampanelli et al.ll2007at [Gonzalez et al.l 
l2007bl ). Kicks this large exceed the escape velocities of 
even the most m assive galaxies, and w ould thus eject SBHs 
from their hosts ([Merritt et al.ll2004l) . Frequent SBH ejec- 
tions would seem to contradict the tightness of the ob- 
served correlations between SBHs and their host galax- 
ies. Kicks would pose an even greater problem at high 
redshifts, where typical galactic escape velocities decrease 
while recoils remain a fixed fraction of the speed of light. 

How might we avoid black-hole mergers that lead to 
large kicks? To answer this question, we must take a closer 
look at how the predicted recoils depend on the dimension- 
less spins Xi = Sj/m? and mass ratio q = 7712/7711 < 1 of 
the merging black holes. Reliable NR simulations have 



been performed for q > 0.1 and \Xi\ < 0.9; in this 
range the recoils are well described by the fitting formula 
(|Campanelli et al.ll2007al) 

v(q,Xi,X2) = v rn&i + v±(cos£ §1 + sin £ e 2 ) + V\\e z (1) 
where 



Br]) 



(2) 



H V 2 A^e z , (3) 

U|| =#/7 2 cos(e-e )|A- L | . (4) 
Here (ei, §2, e z ) are an orthonormal basis with e z parallel 
to the orbital angular momentum L, 77 = q/(l + q) 2 < 1/4 
is the symmetric mass ratio, and A"' 1 " are the components 
of 

A = qX2 - Xl (5) 

parallel and perpendicular to L. O is the angle between 
A. ± and the separation r of the two black holes "at 
merger" . NR simulations indicate that the best-fit val- 
ues for the coefficients appearing i n the above formula ar e 
A = 1.2 x 10 4 km/s, B = -0 93 dGonzalez et al.ll2007a[) . 
H = (6.9 ± 0.5) x 10 3 km/s (lLousto fc Zlochowerl 12008 ). 
and K = (6.0±0.1) x 10 4 km/s (|Campanelli et alJl2007bQ . 
The angle £ ~ 145° for a wide range o f quasi-circular con- 
figurations (jLousto k, Zlochowerll2008lh while 8p depends 
on th e mass ratio q but not the spins ( Lousto fc Zlochowerl 
2009(). The large value of K implies that the merger of 
equal-mass black holes with maximal spins pointed in op- 
posite directions in the orbital plane generates a recoil of 
if/16 = 3, 750 km/s. 

The most obvious way to avoid these large kicks is to 
require black holes to be non-spinning: 

X l = -> A = -> v ± = v\\ = . (6) 
In this case, v m is maximized at t he modest value 175 ± 1 1 
km/s for a mass ratio q ~ 0.36 ()Gonzalez et al.ll2007al ). 
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However, theory shows that non-spinning black holes can 
be spun up to the Kerr limit \x\ = 1 by steady ac- 
creti on after incre asing their mass by only a factor of 
V6 (|Bardeenlll970| ). Observations of Fe Kg fluorescence 
can b e used to measure black-hole spins ([Reynolds et al.l 
1999), and indicate that real SBHs can approach this 
limit: for the SBH in the Seyfert galaxy MCG-06-30-15 
the measured spin is |y| = 0.989tn nno at 90% c onfidence 
(|Brenneman fe Revnoldsll2006t IBerti et al.ll2009ft . 

If black holes are highly spinning, the recoil can be re- 
duced by aligning the spins and thus A with the orbital 
angular momentum: 



X 4 ||L^A^=0^ V ||=0. (7) 

This spin configuration leads to smaller kicks, since the 
coefficient H is almost an order of magnitude less than 
K. The merger of equal-mass, maximally spinning black 
holes with one spin aligned with L and the other anti- 
aligned generates a recoil of H/1Q = 431 km/s. Gaseous 
accretion disks are needed to provide dynamical friction 
to allow SBHs separ ated by r ~ 1 pc to merge in less 
than a Hubble time (Bcgclma n et al.lll98 0l). These same 
accretion disks can exert torques on the SBHs which align 
their spins and orbital angular momentum with that of 
the disk, thus producing the de sired aligned spin config- 
uration ([Bogdanovic et aI1l2007t ). However, the effective- 
ness of this alignment mechanism depends on the highly 
uncertain na ture of the a ccreti on flow near the merging 
black holes. iDotti et ail ([2010) find a residual misalign- 
ment of 10° (30°) between the black- hole spins and ac- 
cretion disk depending on whether the disk is cold (hot). 
This misalignment could be even greater in a gas-poor 
merger or one in which accretion onto the SBHs proceeds 
through a series of small-scale, randomly orie nted events 
(jKing fe Pringlell2007l: IBerti fe Volonterill2008D . 

In this Letter, we present a new mechanism to re- 
duce gravitational recoils by aligni ng black-ho le spin s 
with each other prior to merger. iBovle et"ail ([2008); 
iBovle fe Kesd en (2008) showed that the symmetries of bi- 
nary black-hole systems imply that recoils are only gen- 
erated by a weighted difference of the two spins. This 
general result can be seen to hold for the fitting formula of 
ICampanelli et al.l (|2007aft by noting that a weighted differ- 
ence of spins appears in the numerator of A in equation 
([S]). Spin alignment is a consequence of relativistic spin 
precession as the black holes inspiral due to the loss of 
energy and angular momentum to gravitational radiation. 
We begin calculating the inspiral at an initial separation 
ri = 500i?5 w here spin alignmen t begins for comparable- 
mass binaries (jSchnittman Il200"4l ). and end at a final sep- 
aration rf — 5Rs near where NR simulations typically 
begin. Here Rs = 2GM/c 2 is the Schwarzschild radius of 
a non-spinning black hole of mass M . Relativists use units 
in which G = c = 1, allowing them to measure distance 
and time in units of M , where M = mi + mi is the sum 
of the masses of the merging black holes. We shall do this 
for the rest of the paper. The spin alignment discussed in 
this Letter occurs for both gas-rich and gas-poor mergers, 
as gravitational radiation (GR) dominates the dynamics 
even in the presence of gas at binary separations less than 

™~ 3000m ' 1/4 (im^Ft) ' • < 8 > 



where M is the rate of gas infall (jBegelman et al.lll980h . 

We briefly describe the relativistic dynamics leading 
to spin alignment in § [5] readers interested in further 
details can find them in our longer paper on how spin 
alignment affects th e distributions of black- hole final spins 
( Ke sden et al.l2010t ) . The most notable effect of spin align- 
ment is to suppress the recoil velocity when the spin of the 
larger black hole is initially partially aligned with L. The 
magnitude of this suppression for distributions with differ- 
ent mass ratios and initial spins is presented in § |3J Some 
concluding remarks are provided in § 2J 

2. SPIN ALIGNMENT 

To understand why black-hole spins align, we must first 
describe how they precess at separations r < tqr where 
the inspiral is determined predominantly by the loss of 
energy and angular momentum to gravitational radiation. 
When the orbital speed is much less than the speed of light, 
the motion of the black holes, the precession of their spins, 
and the emission of gravitational radiation can all be cal- 
culated in the post-Newtonian (PN) limit. We us e the PN 
equati ons for precessing binaries first derived in iKidder I 
(|1995|) . supplemented wi th the q u adrup ole-monopole in- 
teraction considered in iRacine I (|2008|) . Gravitational 
radia tion circularizes eccentric orbits ([Peters fe Mathews! 
fl96l . so we restrict our attention to quasi-circular orbits 
of slowly decreasing radius. In the PN limit r ^> M, the 
orbital period i or b oc r 3 / 2 (Kepler's third law), the spin 
precession period t p ~ | S~2^ | 1 cx r 5 / 2 , and the radiation 
reaction time t r ~ r/f oc r 4 . These scalings imply that 
allowing several simplifications. Since 



t, 



orb 



< t p < t r 



(9) 



torb *C t p , the spins S^ process according to 
— " — — — S~2 ? " x Si , 

dt 

where tti are the orbit-averaged spin precession frequen- 
cies. These depend on the mass ratio q, the orbital angular 
momentum L , and the spins Si and S2 as shown in equa- 
tion (2.2) of iKesden et al.l (|2010| ). Since t v <c t r , the total 
angular momentum J = L + Si - 
are constant on the timescale t p 
direction of L evolves as 

dt _ 1 /rfSi 
~dt ~ ]L| \ dT ' 
The magnitude |L| does decrease on the longer timescale t r 
as the o rbital frequen c y uj in creases according to equation 



S2 and magnitude |L| 
This implies that the 



dS 2 
dt 



(10) 



(2 6) o f Kesde n et al 



I2OIOD . 

ISchnittman I ( 2004 ) discovered that if black holes in- 
spiral as described above, spin precession and radiation 
reaction will act to align Si and S2 with each other if 
61 < 62, where 9\ (62) is the angles between L and Si 
(S 2 ). Conversely, if 9\ > 9 2 , Si and S2 will become anti- 
aligned with each other. This alignment is strongest for 
mass ratios q near unity, though it vanishes for precisely 
equal masses as there is no distinction between Si and S2 
in that case. We show the magnitude of this alignment 
for maximally spinning black holes (|x,| = 1) in Fig. [T] 
and for black holes with \Xi\ = 0.5 in Fig. O The upper 
panels show black holes with the nearly equal mass ratio 
q = 9/11, while the lower panels show the smaller mass ra- 
tio q = 1/3. The black curves show that black holes with 
isotropic spin distributions at = 1000M (flat distribu- 
tions in cos 9 1 , cos 6*2 , A0 = 4>2 — </>i , and cos $12 , the angle 
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between Si and S2) maintain these isotropic distributions 
as they inspiral to rt = lOAf, consist ent with previous 
studies of processing spin distributions rtBogdanovic et al.l 
I2007t iHerrmann et aljlioiol: lLousto et al.ll2009D . Isotropic 
spin distributions will have 6\ < 62 just as often as 6\ > 6 2 , 
implying that just as many spins will become aligned as 
anti-aligned during the inspiral. 

These isotropic spin distributions at r; = 1000M are 
only expected for the most gas-poor mergers; in the pres- 
ence of gas, accretion torques will partially align the spins 
and orbital angular m omentum with that of the disk 
(|Bogd anovic et al.ll2007h . We consider a scenario in which 
the spin Si of the more massive black hole is partially 
aligned with L while the other spin S2 remains isotrop- 
ically distrib uted. This scenari o is c onsistent wi t h tha t 
explored by IChang et all ([20091 ) and iDotti et al.l ([2010D . 
where the more massive black hole or "primary" is at rest 
in the center of the accretion disk while the "secondary" 
migrates inwards. The blue (red) curves in Figs. [T] and [2] 
show the subset of black-hole binaries with the 30% low- 
est (highest) values of 9\ at rj. Since S2 remains isotrop- 
ically distributed at r,-, these subsets are consistent with 
the initially flat distributions of cos 6*12 and A0 shown by 
the horizontal dotted lines. However, the distributions of 
cos 6*12 and A<j> no longer remain flat as the black holes 
inspiral to rj as shown by the solid blue and red curves. 
Those binaries with low (high) values of 9\ at rj have val- 
ues of O12 and A(j> at rj strongly peaked about 0° (180°). 
This alignment is very pronounced for q — 9/11 both for 
\Xi\ = 1 and \Xi\ = 0.5, but it is much less significant for 
9 = 1/3. 

3. KICK SUPPRESSION 

How does this alignment of Si and S2 during the in- 
spiral affect the subsequent recoils? We show histograms 
of expected recoil distributions for maximally spinning 
(\Xi\ — 1) mergers with mass ratios q — 9/11 and 1/3 
in Fig. [31 and for mergers with these same mass ratios 
and spin magnitudes \Xi\ — 0.5 in Fig.HJ The upper pan- 
els of these figures show the recoil distributions for the 
same binaries whose spin alignment was shown in Figs. [1] 
and [21 The dotted curves in Figs. [3] and Fig. 0] show the 
recoils expected if the black holes merged with the same 
spin distribution they had at 7"$ (no spin precession) , while 
the solid curves show the different recoils expected if we 
include the spin alignment that occurs as the black holes 
inspiral from n to rt. Comparing the dotted and solid blue 
curves, we see that spin alignment suppresses the recoils 
expected if accretion torques at separations r > fj have 
partially aligned Si with L (low 9i). The red curves show 
that recoils are boosted if Si is initially anti-aligned with 
L, but we expect this case to be less physically relevant. 

Spin alignment and the subsequent suppression of re- 
coils can be even more effective if initially 6 \ < 3 0° as 
suggested by the simulations of IDotti et al.l ([2010) . In 
the middle panels of Figs. [3] and Fig. we show with 
purple, blue, and green curves the expected recoils when 
6\ = 10°, 20° and 30° at ri while S2 remains isotropi- 
cally distributed. The dotted curves show that even with 
this initial alignment between Si and L, the distribution 
has tails that extend to very large recoils. In Table [TJ we 
show the velocities v$o (wgo) at which the cumulative dis- 



tribution function of the recoils rises above 50% (90%) for 
various mass ratios and spins. The first two entries in the 
first column of this table indicate that even if 6\ = 10°, 
half the recoils are greater than 864 km/s and 10% are 
greater than 1,587 km/s for q = 9/11 and \Xi\ = 1- How- 
ever, if this spin distribution is allowed to precess as the 
black holes inspiral from = 1000M to rj = 10M, the 
resulting recoils are dramatically suppressed. The second 
column shows that for the same mass ratio and spins, 50% 
of the recoils are below 273 km/s and 90% are below 611 
km/s. This reduction is a big effect: most of the black holes 
produced will now remain bound to normal-sized galaxies. 
Notice that the reduction in recoil velocities remains signif- 
icant even for relatively moderate spins {\Xi\ — 0.5). The 
precession-induced recoil suppression is less pronounced 
for lower spins because the spin alignment is reduced and 
the recoil v m due to the mass asymmetry becomes com- 
parable to the recoils v±,v\\ from the spin asymmetry. In 
this regime, ejection from the l argest galaxie s seem s un- 
likely anyway: see e.g. Fig. 2 of lMerritt et al.l (|2004l ). For 
completeness, we show how kicks are enhanced for high 
initial values of 6\ in the bottom panels of Figs. [3] and 
Fig. 21 but don't expect such distributions in astrophysi- 
cal mergers. 

4. DISCUSSION 

In this Letter, we have shown that as black holes inspiral, 
spin precession aligns their spins with each other for the 
spin distributions expected in astrophysical mergers. This 
spin alignment drastically reduces the recoils expected for 
the black holes produced in the binary mergers. Spin align- 
ment is most effective for the highly spinning, comparable- 
mass mergers that are predicted to yield the largest re- 
coils (up to = 3,750 km/s ac cording to the fitting 
formula of lCampanelli et al.l (|2007a[ )). Aligni ng the black- 
hole spins generically suppr esses the recoils (jBovle et al.l 
I2008t iBovle h Kesdenl feOQSI: we found a sim ilar suppres- 
sion w ith the alternative fitting formula of iBaker et alj 
(2008). This spin alignment is a purely relativistic effect 
that will occur for all black-hole mergers, as gravitational 
radiation will always dominate the inspiral for separations 
T < r GR ~ 3000M. As long as torques at r > ran align 
Si and L such that 9\ < 30°, spin alignment during the 
final inspiral will nearly eliminate the v > 1,000 km/s 
recoils that are so difficult to reconcile with galaxies keep- 
ing their supermassive black holes. While there is still 
great uncertainty about how merging black holes interact 
with surrounding gas, the PN spin precession discussed in 
this paper is inevitable and results from well established 
physics. We therefore believe that spin alignment must 
be accounted for in future population studies of merging 
black holes. 

We are grateful to Vitor Cardoso for helping to 
test our numerical implementation of the PN evolu- 
tion. We would also like to thank Enrico Barausse, 
Manuela Campanelli, Yanbei Chen, Pablo Laguna, Car- 
los Lousto, Samaya Nissanke, Evan Ochsner, Sterl Phin- 
ney, Etienne Racine, Luciano Rezzolla and Manuel Tiglio 
for useful discussions. This work was supported by 
grants from the Sherman Fairchild Foundation to Cal- 
tech, by NSF grants No. PHY-0601459 (PI: Thorne) 



4 



Kesden, Sperhake, & Berti 



and PHY-090003 (TeraGrid) and by FCT - Portugal NNX07AH06G (PI: Phinney). E.B.'s research was sup- 
through projects PTDC/FIS/098962/2008, PTDC/CTE- ported by NSF grant PHY-0900735. U.S. acknowledges 
AST/098034/2008 and PTDC/FIS/098032/2008. M.K. support from NSF grant PHY-0652995. 
acknowledges support from NASA BEFS grant 



REFERENCES 



Baker, J. G., Centrella, J., Choi, D. I., Koppitz, M., & van Meter, J. 

2005, Phys. Rev. Lett., 96, 111102 
Baker, J. C, Boggs, W. D., Centrella, J., Kelly, B. J., McWilliams, 

S. T., Miller, M. C, &; van Meter, J. R. 2008, ApJ, 682, L29 
Bardeen, J. M. 1970, Nature, 226, 64 

Begelman, M. C, Blandford, R. D. and Rees, M. J. 1980, Nature, 
287, 307 

Berti, E., Cardoso, V., & Starinets, A. O. 2009, Classical and 

Quantum Gravity, 26, 163001 
Berti, E., & Volonteri, M. 2008, ApJ, 684, 822 

Bogdanovic, T., Reynolds, C. S., & Miller, M. C. 2007, ApJ, 661, 
L147 

Boyle, L., Kesden, M., & Nissanke, S. 2008, Phys. Rev. Lett., 100, 
151101 

Boyle, L. & Kesden, M. 2008, Phys. Rev. D, 78, 024017 
Brcnneman, L. W. and Reynolds, C. S. 2006, ApJ, 652, 1028 
Campanelli, M., Lousto, C. O., Marronetti, P., & Zlochower, Y. 2006, 

Phys. Rev. Lett., 96, 111101 
Campanelli, M., Lousto, C. O., Zlochower, Y., & Merritt, D. 2007a, 

ApJ, 659, L5 
. 2007b, Phys. Rev. Lett., 98, 231102 

Chang, P., Strub be, L. E., Menou, K., & Quataert, E. 2009, MNRAS, 

submitted ( arXiv:0906.0825 ) 
Dotti, M., Volonteri, M., Pcrcgo, A., Colpi, M., Ruszkowski, M. & 

Haardt, F. 2010, MNRAS, 402, 682 



Ferrarese, L. &; Merritt, D. 200, ApJ, 539, L9 

Gonzalez, J. A., Sperhake, U., Bruegmann, B., Hannam, M. & Husa, 

S. 2007a, Phys. Rev. Lett., 98, 091101 
Gonzalez, J. A., Hannam, M. D., Sperhake, U., Bruegmann, B., & 

Husa, S. 2007b, Phys. Rev. Lett., 98, 231101 
Herrmann, F., Silberholz, J., Bellone, M., Guerberoff, G., & Tiglio, 

M. 2010, Classical Quantum Gravity, 27, 032001 
Kesden, M., Sperhake, U., & Berti, E. 2010, Phys. Rev. D, submitted 

(arXiv: 1002. 2643) 
Kidder, L. E. 1995, Phys. Rev. D, 821 
King, A. R. & Pringle, J. E. 2007, MNRAS, 377, L25 
Lousto, C. O., & Zlochower, Y. 2008, Phys. Rev. D, 77, 044028 

- -. 2009, Phys. Rev. D, 79, 064018 
Lousto, C. O., Naka no, H., Z lochower, Y., & Campanelli, M. 2009, 

preprint ( arXiv:0910.3197 ) 
Magorrian, J. et al. 1998, AJ, 115, 2285 

Merritt, D., Milosavljevic, M., Favata, M., Hughes, S. A., & Holz, 

D. E. 2004, ApJ, 607, L9 
Peters, P. C, & Mathews, J. 1963, Phys. Rev., 131, 435 
Pretorius, F. 2005, Phys. Rev. Lett., 95, 121101 
Racine, E. 2008, Phys. Rev. D, 78, 044021 

Reynolds, C. S., Young, A. J., Begelman, M. C, & Fabian, A. C. 

1999, ApJ, 514, 164 
Schnittman, J. D. 2004, Phys. Rev. D, 70, 124020 
Tremaine, S. et al. 2002, ApJ, 574, 740 



Relativistic Recoil Suppression 



5 




-1 -0.75 -0.5 -0.25 0.25 0.5 0.75 1 -180 -135 -90 -45 45 90 135 180 

COS 6 12 A0 



Fig. 1. — Histograms of cos $12 (left) and A(f> (right) for distributions of black-hole spins, where 612 is the angle between Si and S2 and 
A0 is the angle between the projection of these spins onto the orbital plane. The horizontal dotted lines give the distributions at the initial 
separation n = 1000M, while the solid curves show the distributions after those black holes have inspiraled to the final separation r/ = 10M. 
The black holes are maximally spinning (|Xil = 1) an d have mass ratios q = 9/11 (upper panels) or q = 1/3 (lower panels). The black curves 
correspond to initially isotropic distributions of the spin Si of the more massive black hole, while the blue (red) curves show subsets of this 
distribution with the 30% lowest (highest) initial values of 81, the angle between Si and L. 




Fig. 2. — Histograms of the same quantities shown in Fig.[T] except the black holes now have initial spin magnitudes \xi\ = 0.5. 
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Fig. 3. — Histograms of the recoil velocity v for maximally spinning (\Xi\ = 1) black hole mergers with a mass ratio q = 9/11 (left) and 
q = 1/3 (right). The dotted curves show the recoils expected if the black holes merge with the spin distributions specified at the initital 
separation = 1000M. The solid curves show the recoils expected if the spins precess as described in §[2] as they inspiral from ri to the 
final separation i*/ = lOAf prior to merger. The black curves in the top panels correspond to isotropic spin distributions for both black holes, 
while the blue (red) curves show the subsets of the spin configurations with the 30% lowest (highest) values of 9\, the angle between the spin 
of the more massive black hole and the orbital angular momentum. The purple, blue, and green curves in the middle panels show the recoil 
velocities expected if initially 8i = 10°, 20°, and 30° respectively and S2 is isotropically distributed. The yellow, orange, and red curves in 
the bottom panels show the recoil velocities expected if initially 81 = 150°, 160°, and 170° respectively and S2 is isotropically distributed. 
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Table 1 



Velocities v S o («9o) in km/s at which the cumulativ e distribution function for black-hole recoils predicted by 
the rit fitting formula (jcampanelli et al"1l2007ah rises above 50% (90%). entries in each row correspond to 
distributions with the same spin magnitudes and mass ratio q = 7712 /mi. the columns indicate the angle 8l 
between the spin of the more massive black hole and the orbital angular momentum at ti = 1000m, and the 
binary separation {ji = 1000a/ or tf = 10m) at which the spins were used to predict the recoil. 
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